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1. Abstract 


This technical report summarizes findings on the 
adaptive control of a transparent Bridgman cr V stal 
furnace. The task of the process controller is to *stab 
a user specified axial temperature profile by controll g 
the temperatures in eight heating tones. The furnace 
controller is built around a computer. Adaptive 
(Proportional Integral Derivative) and Pole Placemen 
control algorithms are applied. The need for adaptive 
controller stems from the fact that the tone dynamics 
changes with respect to time. The controller tested 
extensively on the Lead Bromide crystal growth. Sever 
different temperature profiles and ampoule's translation 

rates are tried. . 

The second objective of this study is to determine 

feasibility of solid liquid interface quantification y 
image processing. The interface is observed by a color video 
camera and the image data file is processed to determine if 
the interface is flat, convex or concave. 


2. Objectives and expected results 


The objective of this experimental study is twofold 

A. Design, build, and validate a computer based temperature 
controller for a transparent Bridgman crystal growth 

furnace. . 

B. Quantify the solid liquid interface of the crystal by 

image processing techniques. 

The transparent furnace consists of eight electnca y 
heated zones which are not insulated from each other. Due to 
the lack of insulation, heat exchange takes place between 
different zones. The process controller is designed in such 
a way that the zone temperature is kept at the desire 
level , despite the heat exchange effects of the neighboring 
This calls for an adaptive controller design. The 


zones • 


adaptability is implemented in a sense that the c ° ntro1 ^ 
tunes itself to changing dynamics of the heating zone, 
change in dynamics is also caused by the ampoule « hlch 
vertically within the furnace. The controller is fo ™* d b * 
current dynamic model of the heating zone. This model take 
into account the effects of heat exchange between the 
heating zones. This undesirable heat exchange can e 

considered as the process disturbance. The task o e 

adaptive controller is then to keep the zone temperature at 
the reference level despite the effects of changes in the 
zone dynamics and the disturbances. 

The second objective of this study is to quantify the 
shape of the solid liquid interface while the crystal is 
growing. Quantification is understood in a sense that the 
interface position is determined with respect to some fixed 
reference coordinate frame. This information can be used for 
the following purpose. 

a - The effect of process parameters on the crystal growth 
can be studied. In a typical furnace, these parameters are 
the axial furnace temperature profile and the translati 
motion of the ampoule. 

b - The interface profile can be used as a feedback control 
signal. The desired interface shape can be compared with the 
actual interface profile and the resulting error sig 
employed as input to process controller. 

c - The position of the interface, determined through image 
processing, permits us to determine the crystal growth rate. 

Quantification of the interface is performed by the 
color image processing techniques. 


3. Introduction 

One common way to grow crystal is to subject the molten 
material to a temperature gradient. This requires keeping a 
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• for a large diameter Gallium Arsenide crystal 
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4. crystal growth control 
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- Temperature measurement or on 

established in comparison wrth for example P 

U " e the te ” PeratUre 
control problem. 

4.r Single input single output adaptive temperature control 

This mode of control starts ^Tnput 

output model of the heating zone. and y 

r.^tr;rr.r.r r— - - - 

the following form: 

y(t)= .X Y (t-1) +a 2 y(t- 2 , + a n y(t-«> + 

bl u(t-l) +b 2 u(t-2)+b n u(t-n)+v(t) 

V(t)= ox V(t-l)+c 2 v(t-2) + c n v(t-n)+e(t) 

. _ for the heat exchange 
where the term v(t) accounts white 

interactions between the neighboring e 'j lmates o£ 
noise. For adaptive control ° obtained by the Least 

process parameters s Technique. At each 

squares or the Genera ize are subst ituted into 

sampling time, these estima values for the 

the control law as if they are 

unknown parameters ( * i,b i'* nd _° 1> 199 ob] ha ve applied 

[Batur et al 

Single input single “"rameters 

an eight zone tran ^ pa “ n sequen tial Least Squares 

(ai,bi> ar6 d two "different control laws namely adaptive FID 
Technique and two controllers have been applied, 

and adaptive pole placem control laws are 

Detailed derivations of these 
given in the Appendix 1. 



4.2 Multi input multi output adaptive temperature control 

Figure 2a shows structure of the 
multivariable controller. Process outputs are the measured 
zone temperatures. Process inputs are the heating powers and 
the translational velocity of the ampoule with respect to 
the furnace. Due to the multizone configuration , there are 
inevitable interactions between the neighboring heating 
zones. The input output dynamics of the system can be 
represented by the following first order multivariable 

model . 

Yi(t)= an yi(t-i)+a 12 y 2 (t-D + • • - +a i, 8 

bn U 1 (t-l)+b 12 u 2 (t-l)+...+bi,9 U 9 (t-1) + 

On e 1 (t-l)+c 1 2 e 2 (t-1) +. • • +c l , 8 e 8( t_1 ) 

y 2 (t)= a 21 yi(t-l)+a 22 y 2 (t-1) +• . -+a 2 , 8 y8( t " 1 ) + 
b 2 i u 1 (t-l)+b 22 u 2 (t-l)+. • -+b 2 , 9 u 9 (t-l)+ 
c 2 i e 1 (t-l)+c 22 e 2 (t-1) +• • • +c 2 1 8 e 8 (t-l) 

y 8 (t)= a 8 1 1 yi(t-l)+a 8 , 2 y 2 (t-l)+ +a 8,8 Y8( t 1 ) + 

b 8/1 u 1 (t-l)+b 8/2 u 2 (t-l)+. • -+b 9/ 9 u 9 (t-l)+ 
c 8 ,i e 1 (t-l)+c 8>2 e 2 (t-1) + +c 8/8 e 8 (t-1) 

These equations can be expressed in a much more compact form 
by the use of vector matrix representation. Combining al 
input, output and disturbance terms into appropriate 
vectors, the multi input, multi output representation takes 

the following form. 

y( t) = A(z- 1 )y(t) + B(z- 1 )u(t-i)+[i+c ( z- 1 )]e(t) (1) 

where the z" 1 is the delay operator defined as, for example, 
y (t-1) = Z _1 y (t) . 


Input output and disturbance 


vectors are explicitly given as 


u(t) 
y(t) 
e (t) 


[ui ,u 2 
[yi i y 2 

t e l ' e 2 


, . • -u 9 ] 
, . • -Y8] 

i • • • e 8 ] 


T 

T 

T 


where T is the transpose operator and u 9 corresponds to the 
translational speed. The A, B and C of (1) are appropriate 
matrices representing the dynamics of the system an 
noise. For example, the A matrix is 


an z 1 ... a l , 8 z 

A(z -1 ) = a 21 z 1 ... a 2,8 z 

ag 1 1 z -1 ... a 8 , 8 z 

The rationale behind this model are as 

follows: 

a- open loop step response of each indivxdual zone 
essentially follows an exponential curve suggesting first 

order dynamics. # 

b- interaction between the heating zones can be taken in o 

account with the off diagonal terms of the A,B,C matrices. 
Note that this multivariable model is drastically different 
from that of the single input single output. Zone 

interactions are explicitly modelled in the multivariable 
case as the off diagonal terms of A,B, and C matrices. The 
single input single output model on the other hand tries to 
accommodate this interactions by the added noise term to 
model. If there is negligible interaction, then A,B and 
matrices become diagonal as , for example 


A ( z -1 ) = {aiiz _1 } 8 *8 i— 1 , 2 . . 8 . 

If the disturbances are not correlated i.e. C={0} then the 
identification of model parameters A and B can be perfor 
by the sequential Least Squares parameter estimation 


techniques. In the case of correlated disturbances, C-{0}, 
The Maximum Likelihood or Extended Least Squares 
identification procedures must be used [Ljung and 
Soderstrom 1983] , [Goodwin and Sin 1984]. 

Once the feedback control system 
performance has been specified, the control action for each 
zone can be determined. In crystal growth furnaces, due to 
the slow motion of the ampoule with respect to the furnace, 
the dynamics of the furnace changes in time. Therefore A,B, 
and C matrices describing the dynamics are not actually time 
invariant. However these changes are not expected to be 
fast. This is because the relative speed of the ampoule is 
slow (on the average 2 cm/day) . Therefore a self tuning 
multivariable temperature control algorithm can be applied. 
Process parameters A, B, and C of (1) can be estimated 
by an on-line Generalized Least Squares algorithm and the 
following performance indices can be minimized by the 
controller . 

4 . 2.1 Multivariable minimum variance controller 


The minimum variance controller minimizes 

V= ( 1/N) { W X (T d (k) -T (k-l-m) ) +W 2 (u (t) ) } ( 2 ) 

where V is a scalar function of input u(t) , output T(k) , and 
the desired process output T d (k) , and (k) is the time index. 
The functions W1 ( . ) and W2 ( . ) define the emphasis that one 
puts on the output error [ T d (k) -T(k-l) ] and the control 
signal u(t) . Common choices are: 

1. W 2 (u(t)) = Ui 2 (t) +. . .+u 9 2 (t) 

2 

2. W 2 (u (t) ) = [u 1 (t)-u 1 (t-l) ] 2 + . . .+[u 9 (t)-u 9 (t-l) ] 



3. 


W 1 (T d (k)-T(k-l)) = [T dl (t)-T 1 (t-in 2 +--- + tT d 3 ( t)- 


T 8 (t-1>] 2 + 

where the choice (1) considers the heating power and 

penalizes excessive power sent to heaters wher 
performs the same action for changes in heating perns . 
inclusion of W 2 ( . ) into the performance index stems from 
the fact that for a typical crystal growth furnace, 
radiative energy received by the ampoule may be a 

considerable portion of the total energy (convective 
radiation). For example, in one experiment invo ving 
eight zone furnace, it has been observed that the hea ing 

wire surface temperature can approach 800 C. At the 

emissivity of 0.8 , approximately 40% of the radiation 

energy will transmit through the quartz ampoule and rea 

the material. . . ^ 

„ l( .) shown in the third choice makes sure that controller 

tracks the desired temperatures in a mean square sense. 


4.2.2 Pole placement control algorithm 

ignoring disturbances, the process dynamics 
of Eqn.l can be expressed as 


y (t) = P(z 1 )u(t) 

If the controller dynamics is written as 


u(t) = C ( z -1 ) [T d (t)-T(t)] 

then the closed loop system, as far as the error is 
concerned, will be given as 


[T d (t) ~T ( t) ] 


{I- [I+P(z - 1 )C(z x )] 


,-ln-l 


P(z -1 ) ] 


C ( z -1 ) ] } 


r (t) 


where I is an 8*8 unity matrix. In pole placement design the 

controller matrix CliT 1 );^ chosen such that thE err ° r 

recovery follows a user specified pattern. 

The reason behind these different choices 

is that each design emphasizes on certain aspects o£ the 
control. The minimum variance controller provides the 
prediction for the ..predictable., part of the disturbances 
and tries to eliminate the disturbances before their effects 
actually take place. The pole placement design however, 
shapes the dynamics of error recovery. Here the disturbances 

are not explicitly predicted. 

in conclusion, the self-tuning control 

algorithm, , _ .. 

a *i A B and C matrices ot rne 

- experimentally identifies a,b, 

process dynamics in real time, and 

- generates a control action for each zone in such a way 
that the performance index in Eqn. (2) is minimized. 

in this report, only the experimental 

results associated with the single input single output 
adaptive control will be presented. Figure 2b shows the 
structure of the single input single output adaptive 
controller and its associated hardware. Multi input ou pu 
experimental work is currently being performed. 

5. Electro dynamic gradient control (ED6) 

Electro dynamic gradient controller 
translates the temperature gradient, instead of the ampoule, 
in principle, this can be achieved by changing the desire 
temperature in each zone so that the effect of translating 
the ampoule can be duplicated by the changes in the esire 
temperature profile. There are two main requirements for 
the success of this procedure • 

1 When the desired zone temperatures are changed, the 
actual zone temperatures should settle on the new levels as 
soon as possible. Furthermore, there should be no 


temperature overshoot during this transition. Due to 
continuous changes in the desired temperatures, this is no 
longer a temperature regulator problem. In controls 
terminology, this is known as a servo problem where the 
desired levels change with respect to time. 

2. No matter how sophisticated the controller is, the fine 
vertical motion of the critical temperature of the imposed 
temperature gradient is effectively bounded by the number of 
heating zones surrounding the ampoule. Increasing the number 
of zones will permit to move this point in very fine 
increments, therefore minimizing the effects of step 
temperature changes within the ampoule. 

The proposed adaptive control schemes will 
also work under this servo control condition. Due to slow 
translational changes in the imposed temperature profile , 
process dynamics, therefore A,B and C matrices are not 
expected to change fast. The trade-off between the fast 
settling and the low overshoot can be adjusted by simply 
changing weights in the performance indices. 

6. Interface Quantification Problem 

There are essentially four reasons to quantify the 

interface between the solid and the melt: 

1. Given the shape of the interface one can determine 
the concavity or convexity of the solid with respect to the 
liquid. 

2. It is possible to determine the relation between the 
interface growth velocity and the ampoule's translational 

speed . 

3 . The effects of the imposed temperature profile on 
the shape of the interface can be determined. 

4. Most importantly, the shape of the interface can be 
used as a feedback signal by manipulating on the 
translational speed and the axial temperature profile. 



The interface quantification through imaging has been 
attempted before. For example [Bachman and Kirsh 1970] and 
have used video signals to locate the melt-solid interface. 
[Fripp et al 1988] have devised an an X-ray imaging 
technique . This procedure works on the difference in 
densities between the transmitted images of the solid and 
liquid phases. 

Figure 3 shows the hardware features of the image 
processing system. The color video camera looks at the solid 
liquid interface within the ampoule. The location of the 
interface within the image frame does not significantly 
change in time. In fact , from the camera point of view the 
relative speed v of the interface is determined by 

v= speed of ampoule - crystal growth speed 

For slow translation rate, the crystal growth rate matches 
the speed of the ampoule and therefore (v) is small and the 
interface practically remains within the viewing area of the 
camera . 

The image is digitized with a resolution of 8 
bits/color and on a three color (RGB) basis. The task of the 
image processing algorithm is to 

i) locate the interface, 

ii) quantify the degree of concavity or convexity , and 

iii) determine the relative speed of the interface , 

in almost real time so that the interface growth can 
further be manipulated by changing the translational speed 
of the ampoule. 

The procedure to accomplish the interface location can 
be summarized as follows : 

1. obtain the histogram of the area of interest within the 
image , 

2. perform histogram based segmentation, 

3. perform dilation or median based filtering to eliminate 
the noise in the segmented image and, finally 


4 apply edge detection procedures to find the interface 
between the solid and the melt . [Rosenfeld and Kak 1982] , 
[Jain 1989] , [Ekstrom 1984]. 


7 . Experimental results 

The schematic ot the transparent furnace is shown . in 
Figure 4. Each one of the four zones located in the ml e 

have a nominal resistance of approximately 8 ohms. u e 
zones have 5 ohms resistance. Temperatures in each zone are 
measured by two thermocouples, one for the 
control and the other for an independent On-Off contro 
for safety. Extra thermocouples are also placed radl y 
axially around the interface to monitor the tempera ure 
distribution, however these are not used by the temperature 
controller. Temperatures are measured by prog 
transducers, currently a 0-700 c range corresponds to 
10 V span. However the range can be made narrower f 
increased resolution without affecting the span. On the 
computer side , the data acquisition board has 
providing temperature reading resolution of 0.17 C. ea 
actuators are zero crossing , duty-cycle modulated solid 
state switches . Each can deliver 30 A. at 110 V. e a 
the computer side of the instrumentation can also 

in Figure rlgure 5 shows the temperature profiles 

when step changes are introduced in the J/'™ 
temperatures. In this case, the adaptive con 
single input single output pole placement controller. We 
have found that it is intuitively more appealing to use 
pole placement controller. This is because, the desire 
dynamic response of the controller can be easily fixed b 
placing closed system poles at desired locations. In th s 
example the desired closed system poles are located 


z = 0.9 , z 2 = O.o. Table 1 shows the estimates of the 

parameters representing the zone dynamics. 

The overall performance of the controller 

is close to the design specifications. For steady state 
operations (all reference temperatures are constant) , each 
zone temperature has controller error with zero mean and 
0.1 C variance. For the dynamic response, the temperature 
control system follows the changes in the desired 
temperatures according to the performance index. 

Figure 6 shows the main result of the 
interface quantification procedure. The interface growth 
rate, which is a critical process parameter, is 
automatically determined by the image processing technique. 
The flow chart for finding this interface is shown in 
Figure 7 . Currently this procedure takes about 5 seconds on 
an IBM PC AT based computer. This time is slightly longer 
than controller sampling time which is 2 seconds. For the 
existing hardware, one practical solution, to this problem is 
to perform the interface quantification less frequently than 
5 seconds. This approach can be justifiable due to the fact 
that the dynamics of the interface is slower than that of 
the furnace. 

The listing of the temperature control and the image 
quantification program is given in the Appendix 2. 

8. Conclusion 

A Bridgman type transparent crystal growth furnace is 
controlled by two different adaptive controllers. Changes in 
the furnace dynamics and the interaction between the 
neighboring heating zones are estimated by an on line Least 
Squares algorithm. The pole placement and the PID control 
algorithms are extensively tested on the growth of Lead 
Bromide crystals. Both controllers performed according to 
their expected levels of performance. However it is found 
that the pole placement controller is much more flexible to 


implement any desired dosed loop teedbac* dynamics on the 
other hand, the adaptive PID controller only executes the 

Zieqler-Nichols based performance index. 

The second objective of this study is to determine 

the position and the shape of the solid liquid interface by 
image processing. The interface is observed by a color. video 
camera and the image data file is processed to determrne if 
the interface is flat, convex or concave. 
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Figure 2a. Adaptive multivariable temperature control system 
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Figure 2b. Adaptive single input single output control system 
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Figure 3. Image processing hardware 
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Figure 4. Transparent furnace and the control system hardware 






Figure 6a Unprocessed image 


Figure 6c Noise is filtered out 
from image 6b. 




Figure 6b Histogram based 

image segmentation 


Figure 6d Gradient based interface 
detection 


Figure 6 Stages of preliminary interface quantification 
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Figure 7. Flow chart for the determination of interface 









Table 


Zone 

a 

b 

1 

-0.9900 

0.0171 

2 

-0.9905 

0.0151 

3 

-0.9917 

0.0123 

4 

-0.9911 

0.0138 

5 

-0 . 9905 

0.0149 

6 

-0.9918 

0.0120 

7 

-0.9894 

0.0187 

8 

-0.9920 

0.0146 


Dynamics of each zone identified in open loop 
experiment 

Process model is y(t)= -a y(t-l)+b u (t-l)+n(t) 


APPENDIX 1 

A.1 Process Modeling 

For the following process structure 


A 1 


yU) = ay(t-l) + 6uU-l) + eU) 


= |y0-l) u 0-1)] 


+ eU) 


A1 


= c}) 0 + e(/) 


Process dynamics 0 can be identified sequentially by, [Goodwin 1984] 


P0-1) l) T 

0 (t) = 9 (/- 1) + [y0) - <}> (t-i) 0 (f— 1)1 

l + 1) PU-l) 4>a— l) 


where 


P (/— l) 4>(*-i) 4> (i-D P 0-D 

PW = P0-1) + 

1 + 4> 0—1) P 0 — 1) <J> 0 — 1 ) 


with the initial choices 0(0), P(0) = a I where a is a large number indicating little 
confidence on the initial choices. 


A. 2 Self-Tuned PID Controller 

Following [Astrom and Hagglund 1988], the controller is given by 

Proportional action: P = K [b r U) — y (0) 


Deriua live action : 


Td dD 
N dt 


+ D = 


-KT 


dy 
d dt 


A4 


dl K 

Integral action: — = — |r(/) — y (/)| 
dt T . 


AS 


A2 


where b, N are appropriate constants, Td is the derivative time constant, K is the 
controller gain and (1/Tj) is the gain associated with the integral action. If the 


sampling time is T then (A3-A5) can be discretized to 

P(k) = K[br(k)-y(k)] 


A6 


D(k) = 


7\ + NT 


D(k-l) 


kt a n 

a 

T, + NT 

a 


) 


[y(k)-y(k- 1)1 


A7 


I (k) = / (*-!) + 


KT T 

[r(k) -y(k) } + — [u(*) - i>(*)l 

T. T t 

l t 


A8 


where the last term in (A8) avoids the integration wind-up by subtracting the 
controller output determined by the controller, i.e. 

= p{k) + D(k) + I(k) 


A9 


from the actual control output 

u(k) = Saturation ( v(k )) 


A10 


which goes to the process in the case of saturation. In (A8), T t is a user specified 
parameter which effectively determines the wind-up recovery. 

Once the process parameters (a, b) are identified, the process gain (k) and the 
time constant (x) can be determined from 


y(k) 
u ( k ) 


bz _1 ru“') , 

<=> — = / 

1 — a 2 ~ * u (z ) 


1 -e 


-sT 


k 


s is+ 1 


= k 


1 -e- TU z~ l 


All 


which yields 


b = k(\ -e~ Tlt ) , 


a — e 


-T!\ 


A12 


or 


A3 


k = , i= —T/(na 

1 + a 


A13 


Following the Ziegler and Nichols tuning rules, the controller parameters can be 
determined by 


K = 1 . 2 — ; T.-2T ; 7 *. = 772 
kT ‘ d 

where it is assumed that the process delay time is equal to the sampling period. 


A.3 Self-Tuned Pole Placement Controller 

The pole placement controller is assumed to be in the following form 


u ( k ) = — 

) 




A14 


where R (z- 1 ), S (z-i) and T (z-i) are polynomials to be determined. Writing the 
process dynamics from (A1) as 


yjft _ bz~ l = BU *) A15 

u (k) 1 — az -1 

The dynamics of the feedback control system is given by 


y(k) = 


BT 

AR + BS 


r(k) 


A16 


where for convenience the argument (z-i) in the polynomials has been dropped. If 
the desired characteristic polynomial for the feedback control system is given as P, 
then from (A16), it follows that 


AR + BS = FP 


A17 


A4 


where F is an observer polynomial. To guarantee the integral effect, the polynomial 
R is designed to have an integrator, i.e. 


Under this condition (A17) becomes 


(1— az 1 ) (1 — 2 l )(l+qz ] ) + bz 1 (s Q + Sj ^ 1 ) = (1 +fz 'HI+Pj* 1 + p^z 2 ) /\ig 

Equating coefficients associated with the equal powers of z gives, 

~fp 2 f + p } + a + 1 - (fp 2 /a) f Pl + p 2 - a - {-1 - a)(fp 2 /a) 

9= __. So= ; s 1= - 

under this solution and from (A1 6) and (A 17) the closed loop feedback control 
system can be represented by 


Bir'inr 1 ) 

y(k) = — r(fe) 

F(z- 1 )P(2 _1 ) 


A21 


if 


then (A21) becomes 


T(z~ x ) = F(z~' 


/Ml] 

in D 


B(z _1 )P(l) 

y(k) = ; r(k) 

P(z~ l )B( 1) 


A22 


which yields 


- _ A23 

y = r 


APPENDIX 2. Ccfntrol and image processing software 


{iBclude<stdio.h> 

{include<ctype.b> 

|include<dos.b> 

{include<graphic3.h> 

|include(natb.b> 

{include(conio.h) 

{include<tardev.h> 

{include ’coi.h" 

{define live 1 

{define width 0.75 /* anpule width in inches */ 


int Day0.PlotI-301; 


uinfl 

{ 


unsigned long 
int 
char 
int 
void 
void 
extern 
void 
void 
void 
void 
char 
void 
void 
void 
void 


TinelntoTest (unsigned long StartTine); 
sure (ch«r[201 ) ; 

GetRef (); 

KTD(int Channel); 

d2a(int dio.int channel, int base); 

ZeroOut(void) ; 

void Drawlxis(void); 

PlotLine (float Y(8], float U[8], float Phi 181 [2], int HPORT); 
cl (void) ; 

■enulint, int, int, int) ; 
clearshov(void); 

getref ( long int SlTine.long int RTine); 
ccob( ) ; 

getline (unsigned char b[ 5121 . unsigned char *no); 

heunkel (unsigned char searchln[512), unsigned char ‘no, unsigned char ref In 15121 ) 
drawln(unsigned char v(256], unsigned char len); 


struct tiie now; 
struct date dnow; 

float c [8] 121 = 11-5.4395, 0.9867), (1.9618, 1.01), (-6. 7497, 0.99171, (8. 9301 .0.9869), 

(-1 6935, 0.9979), (0. 6859.0.9744). (-1.7250. 1.0032), (2. 2191, 0.9964)), 

Phi 181 12 ) , Theta 18)121, P18) 12) 12), PxPhi (21, Yl81,Ul81,r 181, ref 18), ry 181, 
den, res, Unax, delh, input, 

sO, si, rlu,rbar,polel,pole2. alpha, d.incden.UDelay 181, YDelay 181, 

terp; 

float tspeed.anpulescale; 

char Status, cad, filen|121 , constringllOO) ; 

unsigned long StartTine, RTine, SITine.DataTine; 

unsigned char i.j.k, Done, HPORT. Plotl.ys.xs, Simulate, Save. ReDrav.shovo, 

datai, identify, h.n.s.hs.shovid, interface, tenpc, 
search In 151 2) , lnlen.ref In 1512) , ret In len, ref of fset, trans; 
int novi, dwell, tdist , delh2; 


FILE ‘f; 



Sinulate = 0; 

ReDraw = 0; 

show * 0; 

should = 0; 

identify * 1; 

directvideo ■ 0; 

interface = 0; 

trans = 0; 

textcolor(3); 

Uwx * 600.; 

NPORT * 7; 

Save = 0; 

Done * 0; 
delb -5.; 
delh2 - 100‘delh; 
data! ■ 1.; 

Plotl - 301; 

RTik = 0; 

SITiie » 0; 

polel * -1.5; 
pole2 • 0.56; 
d - l*polel+pole2; 
alpha = -0.5; 
incden* 0.00001; 

Graphlnit(-l); 

SetPageHode(4); 

if (live) SetBLiveModeO ; 

else 6etPic(’furn256.tga\ 0,0,0); 

for(k-0;k<-HP0RT;kH) 

{ for(i=0; i<=l ; i++) 

for(j-0;i<*l; jt+) 

{ Phi Ik] li) -1.; 

Theta [k] 10] =-0.9; 

Theta Ik] (11 =0.06; 

Plklliim *0.; 
if(i-i) Plkllilli) - 1000.; 

1 

UDelaylkl = 0.; 

YDelaylkl = 0.; 
rlkl =1.; 

Ylk] ■ 1.; 

U(kl =0.; 

) 

/* get settings */ 

DravkxisO; 

nenulSave.shovo.showid, interface) ; 


/* zero outputs */ 

if(!Siiiulate) ZeroOut ( ) ; 



do 


StartTi»e=TiiieIiitoTest (0) ; 
SITinet-l; 


/* READING TINE V 

/‘ REFERENCE TEMP */ 


Status-getref (SIT]Be*delh,RTiBe) ; 
svitch(Status) ( 

case T: 

for{j*0;]<*NP0RT;jH) reI[ 3 l=ryUl+(rljl-ry[]])*SITiBe/RTi«e*delh; 
gotoxy(35,9); 

h‘(RTme-SITiBe*delh)/3600, ; 

B c (RTi*e-SIT)*e t delh-h‘3600. 1/60. ; 

s-(RTiBe-SITiBe*delh-h*3600.-B*60.) ; 

cprintf {‘REMAINING RAMP TIME -) X02d:X02d:X02d\b.B.s) ; 

break; 

case *S‘:for(H;j<*NPORT;jH) rellj)*r|il;break; ) 


/* READING TEMPERATURES */ 
/* SIMULATED PROCESS */ 

if (Sinutate) 

for(i'0; K'NPORT; iff) Y{il = -0. 9*Phi [ 1 ] (0]+0. 088823*Phi li] 11] ; 


/‘ ACTUAL PROCESS */ 

if(!Siaulate) 

for (i-0;i(-NPORT;iH) 


{ 

Y[i ] = clil[0Hclil(l]*(ATD(i)tl51. 661/6.0667; 
if (YliJ)lOOO. ) 

{ 

gotoxy{20, 15) ; cprintf ( "Tberoocouple Problea 1)'); 
ZeroOutO; 

Done=l; 


/* 


cprintf("\n‘);‘/ 

/* DETERMINE OUTPUT */ 

for(j=0; ]<=NPORT; jt+1 
{ 

rlu * -1. 0*alplia*pole2/{'n)eta l j) lO]+incden) ; 

sO * (alpba+polel+1. 0-Theta (}] 101 -rlu|/(Theta(jlll]+incden) ; 

si • ({alpha‘polel+pole2)-(Tbeta[]l 101-1. I'rlntThetalfl [01 l/(Tbetal]llll+incden); 

rbar * ref 1 i 1 * (1 . ^a lpha) *d/ (Theta l j 1 [1 1 i-iBcden] ; 

Ulil * (1 . 0-rlu) *U( j]+rlu*UDelay l j l-sO*Ylil-sl t YDelayl jl+rbar; 

if (U(]l> c UBax) U[ j 1 = Uaax; 

if (U[)l<01 UU1 = 0; 


/♦ SEND POWER ‘/ 

i f { .' Siaulate && I Done)) 

for(j-0; j<-5; j++l d2a (Ulil, j, 512); 

for{ j=0; j<“l; jaa) d2a(Uljt61.]H.528); > 


/‘ SA.VE DXU '/ 

it (Save TinelntoTest(DataTiiie) >=datai*100) 

( 

DataTue-TuelntoTest(O) ; 

gettinednov); . .. 

fpriDtt(t,‘W2d:X02d:W2d.W2d’.now.ti_hottr.Dow.ti_iin.now.ti_sec,now.t)JunQ); 
for(i*0;i<=NP0RT;iB) f print! (t. * *7.2f\Ylil); 

it ( interface) 

( 

if (live) 

( GrabFraneO; 

SetDispHodeO; 

) 

else 6etPic(*furn256.tg«” . 0.0.0) ; 

heunke 1 (searchln, ilnlen, ref In) ; 
forti-0; i<lalen; i+*5) 

( te»p=(searcbln[256+i 1-ref ln(256+i fre! of f set ] +0 . 1 ‘anpulesca le ; 

fprintf(f,‘M.3f'.tenp); 

1 

if (live) SetBLiveModeO ; 

else GetPict“furn256.tg#".0.0.0) ; 

) 

fprintf(t,*W); 

* /* IDEHTIFY SYSTEM */ 

if (identify tt rt3))l) 

for(k-0;k<*KPORT;kH) 

{ for(i=0;i<*l;iB) 

( PrPbi l i 1=0. ; 

for (j*0; j<*l; j++) 

PxPbi lil+=P(kl til [ jl*Pbi Ikl lil : 

) 

den*0. ; 

for(H;i<“l;iB) 

den Phi [k] lil*PxPbi lil ; 

den 0.0001; 

for(i~0; i<”l; i++) 
for 

Plkllillil -=PxPbi lil*PxPbilil/den; 

res*0.; 

for(i«0;i<*l;iB) 

res b Phi (kllil ‘Theta (k| lil; 

res -Ylkl-res; 

if( fabs(res) > 0 . 0 ) 

for( i-0; i<“l; i++ ) 

Theta Ik 1 1 1 ) +=PxPhi ( ) 1 /den * res ; 
if (ThetalklllXO ’,! Theta Ik] HI >2) Theta Ikl [1H. 06; 
it (Theta Ikl 101 >0 Thetalkl 101<-1) Theta Ikl 101-0.9; 


/* PLOT DkTk */ 

PlotLine(Y,l),Phi,NPORT); 

whi le(Ti«e!ntoTest (StartTine ) <de lh2) 



wh j leiTinelntoTest (StartTiee ) <de lh2) 


,rlHll; 


if (kbhitO) 

{ 


/* KEY PRESSED ? */ 


ced*toupper(getch() ) ; 
switch(ced) 

{ 

case T : if (sure (“CHANGE T£MP’))( 

f or ( j=0; j<=HP0RT; j++i { 
badscan: gotoxy(35, 9) ; 

if(!iovi) scanf(" kc'.Steepc); 

if (j==0) cpriutf (“ENTER ZONE Id NEW TEMP : \jH); 

else cprintf (“ENTER ZONE M MEW TEMP (tt.lf): “.3+1 


eovi'Scanf(" Xf’.Steep) ; 
if(eovi) rlil'teup; 
rym-YUh 
cl { ) : } 

gotoxy(35,9); 

cprintf (“ENTER RAMP TIME (Mrs): 
scanf (“If “ . Steep ) ; RTiee'teep‘3600. ; 

SITiee^O.; 

cl(l;} 

break; 

case “E* : if (sure(“EHT“) ) 

Done-1 ; 
break; 

case T ; i f (sure (“START/CHAN6E TRANSLATION*) ) ( 

trans-l;gotoxy(35.9); 

installcoedrivers (1); I* Set up for CONI: 


*/ 

N 8 1 */ 


ist=9843./eovi;if(tdist<0) tdist‘=-l; 


setbaud (2400); /* 2400 Baud 

cprintf (“ENTER THE TRAMS DISTANCE (+/- ce); *); 

scanf(“W“,Stdist);cl(); 

cprintf (“ENTER THE TRANS SPEED (CM/DAY): “); 

scanf(‘*f\Stspeed);cl(); 

if (tdist<0) (eovi=-l;tdist*»-l;}else eovi'l; 

dwel l=8778./tspeed; 

while (dwell<500)( 

if(eovi<0) (iiovi— ; dwe 11=8778 . /tspeed*-eovi ; ) 

else (BoviB;dwell=8778./tspeed‘BOvi; } 

1 


coejrintf 

CL0AD;\n\r“); 




coi_pnntf 

(“10MODE M INC 

NULL 

NULL 

\n\r“); 

coejrintf 

(“15MATH Q1 


0 

\n\r“); 

coejrintf 

(“20LOOP 

%i 


\n\r“. tdist) ;td 

coejrintf 

(“30LOOP 

%i 


\n\r“,tdist); 

coejrintf 

(“40MOVI 

td 

0 0 

\n\r'.eovi); 

coejrintf 

(“45DISP LNUN 

BLANK BLANK BLANK 

\n\r“); 

coejrintf 

(“50MATH Q1 - 

01 

+ 1 

\n\r“); 

coejrintf 

(“55NRITE 1 

0 

0 Q1 

\n\r“); 

coejrintf 

(“60DVELL 

Id 

MILLISECONDS 

\n\r“,dwell); 

coejrintf 

(“70ENDLOOP 



\n\r“); 

coejrintf 

(“80ENDLOOP 



\n\r“); 

coejrintf 

(“90DONE 



\n\r“); 

coejrintf 

(•♦;•); 





un installcoedrivers {); 


) 



case T : if(!Save AA sureCSTART AQU1SITKT))! 
badfile:cl(); 

cprintf (“ENTER THE FILE NAME: *); 

scanfCks'.filen); 

f-fopen(filen, V); 

if (f=*NULL) goto badfile; 

gettine(Anow); 

getdate(Adnov); 

f print! (f. “File Started *02d/*02d/*02d at *02d:t02d:*02d.*02d\n\dn 


ov.dajDon,dnow.da_day,doow.da_year, 


now.ti_bour,now.ti_»in,now. ti_sec, now, t i_hund); 
Save-1 ; cl ( } ; gotoxy (35, 9) ; 
cpnntf CENTER THE SAMPLE INTERVAL (s): *); 
scant CM', Adatai) ;cl(); 

DataTine-TinelntoTest(O); 
nenu(Save, showo.sbowid. interface); } 

case ‘L* : it ((Save AA (interface AA sureCRECORD INTERFACE PROFILE* I) 

! ! (Save AA interface AA sureCChange SEARCH and/or REFERENCE Line(s)‘)))( 
1 f (1 i ve ) !6rabFrane ( ) ; SetDispMode ( ) ; ) 


if ((interface !! sureCChange SEARCH Line*))! 
clearsbow();gotoxy(l,l); 

cprintf ("Enter an Approxinate SEARCH Line\n\n\r*) ; 
cpnntf (‘ -nark points tron LEFT to RI6HT\n\r*); 
cprintf (‘ -Use SPACE to nark points \n\r*|; 
cprintf (“ -Press RETURN when done*); 
get) ine(searchln.Alnlen) ; 
aiipulescale-width*2.5d/(lnlen+0. ) ; ) 
i f ( ) interface I! sureCChange REFERENCE Line" ) ) { 
clearshow();gotoxy(l.l); 
cprintf ("Enter the REFERENCE Line\n\n\r“); 
cprintf (* -nark points froi LEFT to RI6HT\n\r*); 
cprintf {* -Use SPACE to nark points \n\r*) ; 
cprintf (“ -Press RETURN when done\n\n\r‘); 
cprintf C *Be sure this line is AT LEAST AS\n\r*); 
cprintf (" LONG AS THE SEARCH LINE*); 
getline(refln,Areflnlen); 
refof fset-ref ln[i]-searchln|01 ; > 
cl () ; clearsbowl) ; interface*) ; 
nenu(Save, showo.sbowid, interface) ; } 
if(live) SetBLiveModef); 
break; 

case T : if (Save AA sureCSTOP AQUISITION*)){ 

fclose(f); 

Save-0; 

nenu(Save, showo.sbowid, interface) ; ) 
break; 

case ‘O’ ; if (Ishowo AA sure ( "SHOW OUTPUT*))! 

shovo-1; 

sbowid-0; 

clearshowO; 

nenu(Save, showo.sbowid, interface) ; ) 
break; 

case T ; if(!showid AA sureCSHON IDENTIFIED PARAMETERS*))! 

shovid-1 ; 

showo-0; 

clearshowO; 

nenu (Save .showo.showid.interf ace ) ; } 
break; 


case 


H* : if (showo U sure ( “HIDE OUTPUT*)) ( 

clearshowO; 
showo*0; 

BenutSave, showo.showid, interface) ; ) 
break; 

case 'J* : if (showid l& sure (“HIDE IDENTIFIED PARAMETERS* ) H 

clearshowO; 
showid=0; 

Benu(Save, showo, sbowid, interface) ; ) 
break; 

) 

) 

f or( j“0; j<=MP0RT ; i++) 

{ 

YDelaylil - Ylils 

UDelaylil = Phi [ i ) 111 ; 

Philillll * U(j|; 

PbiljllO) - -YU1: 

) 

if (showo) 

( 

gotoxy(d, 1) ; cprintf (“Zone Ttarg Tact P.0 '); 
for(i*0; K-HPORT; in) 

( 

gotoxy(5, 2+i) ; 

cprintf (*M (6. If *6. If %8.4f *, i+1. r[i] .Yin .Ulil/4025. *100. ) ; 

) 

) 

if (sbowid) 

« 

gotoxy(6, 1) ; cprintf (" MU Bill '); 

for(i=0;i<=NPORT;iH) 

{ gotoxy (5, 2+i ) ; 

cprintf ( * k8.4f k8.4f",Theta[i|(01,Thetali|[l]); 

) 

] /* DELAY 

if (SITiie==3) delh*TineIntoTest(StartTiBe)/100. ; 

) while(!Done); 

if ('Simulate) ZeroOut { ) ; 
if (Save) fclose(f); 
closegraphO; 


unsigned long TineIntoTest(unsigned long StartTine) 

{ 

struct tine now; 
long TifteNov; 
gettine (Anov) ; 

TiieNovUong) now. t i_hour'360000+( long) now. ti_mn‘6000+ (long) now. ti_sec*100+now.ti_hund*l; 

if (TiBeNow(StartTine && StartTiie !-0) return(Ti*eHow , 24*360000-StartTiie) ; 
return (TiieHow-StartTiBe) ; 



char getretdong SITine, long RTine) 

{ if (SITiue < RTine 66 RTine)0 ) 

else return! ‘S* ) : 

) 

int RTDUnt Channel) 

{ 

unsigned int Base.vordl.vordh; 
unsigned char xl.xh; 

Base-768; 

out port (Base+2 .Channel); 
outport(Base.O); 

/* vbile(inport(Base+8)>7)=-l);*/ 
delay(lOO); 

xl-inportb(Base) ;xh-inportb{Base+l) ; 
vordl=xH0.;vordh=xh+0.; 

I* cprintfCtd %i \r\n",wordl,vordh);*/ 
return ( (wordh< < 4 ) + (word 1 > >4 ) ) ; 


void d2a(int dio.int Channel, int Base) 

( 

int xh.xl; 
xh-dio/256.; 
xl-dio-256.*xh; 
outportb(Base+2*Channel , xl ) ; 
outportb(Base+l+2‘Channel,xh) ; 


void ZeroOutO 

{ 

int i; 

for(j-0; j<-5 ; i++) d2a (0. j , 512} ; 
for(j=0; )<-!;)+♦) d2a(0, j+4,528) ; 


void cl() 

{ 

int i; 

gotoxy(35,9); 

for (i-1; i<48; i++) cpnntfl* “); 
gotoxy(35,9); 

} 

void clearshowO 
{ int i.i; 

for (i-l;i<-9;i*t) 

{ gotoxy (1 , 1 ) ; 

for (j*l; j<35; jw) cpnntfC “); 



int sureichar s (121 ) 

{ char cud^'x’; 
cl(); 

gotoxy{35.9); 
cprintffts (y/n) ? *.s); 

while! cud.'-'Y 1 66 end )»■«*) cwHoupper(getchU); 
cl () ; 

return! (ciid== ‘ Y ‘ ) ? 1 ; 0); 

1 


return { ’ U ' 



void getline(unsigned char xl 15121 , unsigned char 'no) 

{ long color 16] . t=0; 
int i.key.xa; 
unsigned char j, x.y; 

H; 

X'50;y=50; 

do{ 

for(i»-l;i<2;iH) (6etPix(icolor t i+U . x+i ,y, Oh PutPix(&t, x+i , y . 0) ; > 
ford— l:i<2;i+-2) {6etPix(&color d*41 . x. y+i . Oh PutPix(&t . x. y+i . 0) ; > 


whiled kk i!-32 && i!-13) 

i=getch( ) ; 


if (i— 0)t 

for(i— 1; i<2; i+4) PutPix(icolor[i+i|.x+i.y.O); 
for(i— 1; i<2; i+-2) PutPix{S«color l i+41 , x, y+i , 0) ; 
key-getch(); 
switch(key) 

{ case 75: x- r 3; break; /‘ left '/ 

case 77 :x+=5; break; /* right */ 
case 72 :y+=3; break; /* up */ 
case 80: y--5 ; break; /* down '/ 


H 


else if (j==0 fl i *=13) (xf [0]=x; xf f 2561=y; j+t; xa=0; } 

else if (i!=13){key=x-xf Ixahif (key<0) key'=-l ; xf [key+xa]=x;xf [256+keytxa]=y; 

for(H;i<key;j++) 

{ xf ( j+xa ] =xf (xa 1 +]*key/ (xf (key+xa ) -xf lxal+0. ) ; 

xf [256+i4xal=xf l256+xal+(xf 1256+key+xal-xf t256+xal)/(xf Ikey+xal-xf [xa]+0. l*(xf [ j+xal-xf (xal ) ; 

PutPix(&t,xf l j+xa 1 . xf 1256+ j+xa 1 , Oh 

) 

xa+-key; } 

}vhile(i!-13); 

*no-xa; 

) 

void drawln(unsigned char v[2561, unsigned char len) 

{ unsigned char i; 
long t-31«5; 
for(i-0;i<len+l; i++) 

PutPix(lt,v(il.vi256dl.01; 

1 

void »enu(int save, int showo.int showid, int identify) 

{ 

int i,x#in*20,yiiin-21; 

char strllllMOH’T -> Change Teeperature References’, 

*T -> Exit’, 

*T -> Acquire Teeperature Data’, 

”C' -> Close Data File’, 

"O' -> Show Outputs’, 

’T -> Hide outputs’, 

*T -> Show Identification Paraneters’, 

"J‘ -> Hide Identification Paraneters’. 

"L 1 -> Change Quant Ref/Search Lines’, 

"L‘ -> Quantify Interface’, 

’T -> Start / Change Translation’); 


setviewport(280,0,639, 100,1) ; 
setcolor(BLUE); 

setfillstyle( LTSLASH_FILL , BLUE ) ; 
bar3d{10,10, 340,90, 0,2); 
setcolor(VHITE); 


for (y=0;y<K; y++) 

for (x=0;X<N; X++) 

d(N*y+x| - x*sina lpha+(H-y+l } ‘cosalpha ; 

for(i-0;i<N‘N; in) 

for(j-i+l; j++) 

if (d(o|)| |<d(o[i) )) 

{ nini = o [ j ] ; 

o[j] * olil; 
oil] * Bini; 

} 

for(i-0; i<K*N; in) 
l sub = 0. ; 

for(H; J< = 1;)H) 

sub +- collo()ll-Bean; 
if (sub)bux :: i--0)( maxi * i ;»ax ■ sub;) 

} 

if (line) { 

for (X-0;X<N;XH) 

for (y-0;y<N;y++) 

if(H*ytx)Baxi) 

( Bini=searcbln[256+k] ; 

searchln|256+k]-ygiy; 

color-0; PutPix(icolor. xg+x.yg+y.O) ; 

color-(long) 255<<16; PutPix(icolor, xgix, ref lnlxg+x+ 256-ref ln|0] 1.0); 

i f (x—ll-1 ) yg-searcb 1 n 1 256+k I -H/2 . ; 

y-H;k++; 

)) 

else] 

for (X'0;X(H;XH) 

for (y-0;y<H;yH) 

( if (d‘y+x(=Baxi) { color=31<<5 ; PutPix(4color. xg+x,yg+y, 0) ; ) 

else (color-0; PutPix(t.color,xg+x,yg*y,0);) 
if(N*y+x>Baxi &£ x—N-1) 

( y»in=yBin+y-H/2. ; 

y-N; 

) 

)) 

) 

) 

void PlotLine(float Y[8] .float U181 .float Pbi (8) 12) . int NPORT) 

I 

int xs.ys, i, j.cl81=(2.3. 4, 15, 13,12,11.14} .color; 

float posl8Hl. 5. 4. 5, 6.5. 8. 5, 10.5. 12. 5. 15. 5. 18. 5); 

static int Plotl , ReDrav.yoldlfl] ; 

if (PlotI>263) /* Erasing TEMP v Tine */ 

( setviewport (56,150,319,323,1); 

clearvievportO; 

Plotl-0; 

setlinestyle(4, 0x8080,1); /* Grid Lines */ 

setco lor(4 ) ; 

for(i-l;i<-7;iH) 

( line(-l, 175-i‘25, 262, 175-1*25); 

yoldlil-0; 

) 

setcolor(15) ; 
setlinestyle(0,2,l); 


setviewport (55,150,319,323.1); 
setlinestyle(4,0xfftf,l); 
color=getcolor(); 

f or (H; IMPORT; )++) /' Plotting TEMP vs TIME */ 

t 

ys - 175-25. /100 . * (Y l j ] +2) +1 ; 
xs - Plotl; 
setcolor(clil); 

1 i ne { xs - 1 . yo 1 d l j 1 . xs . y s 1 ; 
yold(jl=ys; 

) 

ReDraw=0; 
setcolor(color) ; 

tor(i=0;j<-»PORT;i++) if(fabs(Yl]l+Phi U) 101P0.5] ReDrav=l; 

if (ReDraw !! Plotl==0) /* Plotting POSITION vs DISP V 

< 

setviewport (371. 150,639, 314, 0); 
clearviewportO; 

setl inestyle (4 , Oxf f f f ,1) ; 
for(H;j<=NPORT;]H) 

( xs = 35. /100. MYU1+0. 15) ; 

ys = 5tposm*2.54‘32./10.; 
if(i>0) lineto(xs,ys); 
else noveto(xs.ys); 

setlinestyle(4, 0x8080.1); /* Grid */ 

setcolor(lO); 

for(i=l; i<=7; i++) 

line(35‘i, 0,35*i, 165); 

setcolort 15) ; 

) 

PlotI+=l; 


outtextxy(120.yBin-9, “Action Menu - ) ; 
setcolor (WHITE); outtextxy(xmn.ynin+Q*10+50.str 101 ) 

setcolor(RED); outtextxy(XBin.yBin+l 

setcolor(WHlTE) ; 


* 10 + 50 . strlll); 


if (.'save) 
else 

it (,'showo) 
else 

if (Ishovid) 
else 

if (save il identify) 
if (save Ui .'identify) 


outtextxy (xai b , yei n+30, str {21 ) ; 

outtextxy(xBin,ynin+30,strl31); 
outtextxy (xmin, ymn+10, str [41 ) ; 

outtextxy(xnin,ynin+10,str[5] ) ; 
outtextxy(xnin,yBint20. str 1611; 

outtextxy (xnin.yiiin+20. str 17) ) ; 
outtextxy(x*»in.ynin+0, str{81 ) ; 
outtextxy(xmn,yiiin+0,strl91); 

outtextxy(xiiin, yuin+40, strllOl ) ; 


} 

void heunkel (unsigned char searchln 1512) , unsigned char ‘no, unsigned char reflnl512]) 
tdefine N 8 
Idefine H H*H 
( 


unsigned char olMJ . col IM1 . xg.yg.x. y, i . j . bbxi . mini ,ymn. k=0; 
double sinalpha, cosalpha. tannuB. tanden. dWl, pi; 

long color; 

float Bean,suB.iax=0; 

int r.g.b.c.t; 

pi-acosH.); 
yg*searchln{256]-N/2. ; 

for(xg=searchlnlO!;xg<=searchlnl*no| U xg)=searcbln 101 ; xg+=«) 

{ 


*ean*0. ; 

for (y*0;y<H;y++) 

for (X-0;X<N;X++) 

{ GetPixt&color, xg+x.yg+y. 0) ; 

col lM*y+xl ■ ((color»8)&255|; 
Bean += col(H‘y+xl; 

o[H*y+xl * H*y+x; 

) 


nean /* {N*H+1.); 
tannufe - 0.; 
tanden - 0,; 
for (y*0;y<H;y++) 

for (x-0;x<H;x++) 
{ if(x<y) 
else if(x!-y) 
if (x+y>H-l) 
else if(x+y!*H-U 


) 


tannuB += col [N‘y+x] ; 
tannuB — col(N*y+xl; 
tanden += col {H*y+xl ; 
tanden -= col lN*y+x) ; 


if (tanden!=0. ) 

( sinalpha = sin(atan2(tannuB.tanden)-pi/4. ) ; 

cosalpha * cos(atan2(tannuB,tanden)-pi/4. 1 ; 

1 

else (sinalpha * 1.; cosalpha = 0.;) 
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Figure 2. Interface shape control 

(applicable only for transparent furnace and materials) 










